Order parameter of quasi-one-dimensional superconductors: 
symmetry features in quasiparticle density of states and spin susceptibility 
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Recent experiments indicate that the Bechgaard Salts (TMTSF)2C104 and (TMTSF)2PF6 may 
be unconventional triplet superconductors. The quasiparticle density of states and the uniform 
spin susceptibility tensor are computed at low temperatures for order parameter symmetries, as an 
attempt to narrow the number of possibilities based on current experimental evidence. 



In early theoretical investigations of quasi-one- 
dimensional superconductors (of the Bechgaard family 
(TMTSF)2C104 and (TMTSF)2PF6) in the presence 
of an external magnetic field, the possibility of triplet 
behavior in these systems was suggested [Q-^- The 
present experimental evidence indicates that these quasi- 
one-dimensional superconductors are indeed unconven- 
tional 0-01 ■ For instance, Lee et. al. measured 
the magnetic field versus temperature phase diagram for 
(TMTSF)2PF6 under pressure of 6 kbar. They found 
that the upper critical fields along the usual a, b', and c 
direction are highly anisotropic. Furthermore, the Pauli 
paramagnetic limit is exceeded by a factor of 3 or 4. In 
addition, Lee et. al. |0] found that there is no Knight 
shift in (TMTSF)2PF6 for fields H || b' at pressures P « 
6 kbar. This suggests the existence of a triplet supercon- 
ducting phase in this system. The temperature versus 
magnetic field phase diagram of (TMTSF)2C104 at am- 
bient pressure was also measured by Lee et. al. 
for H II 6'. These measurements seem to indicate that 
the Pauli paramagnetic limit is exceeded in this com- 
pound, which suggests triplet superconductivity. Fur- 
thermore, Belin and Behnia |^ reported measurements 
of the thermal conductivity in the superconducting state 
of (TMTSF)2C104, indicating their data is inconsistent 
with the existence of gap nodes at the Fermi surface. 
These experimental results combined suggest the exis- 
tence of a fully gapped triplet superconducting state in 
(TMTSF)2C104. Ahhough the similar crystal structure 
of these systems suggests, from a simple group theoreti- 
cal point of view, that the origin of the pair interaction 
is the same, the role of the proximity to an SDW phase 
in (TMTSF)2PF6, for instance, needs to be investigated 
both theoretically and experimentally. Furthermore, the 
order parameter symmetry also needs to be investigated 
via phase sensitive experiments like those performed in 
cuprate oxides P,p^. 

In this paper we are concerned with the symmetry of 
the order parameter of a triplet quasi-one-dimensional 
superconductor at zero magnetic field and we have three 
main goals. First, we perform a group theoretical analy- 
sis of the possible symmetries of the order parameter for 
an orthorhombic quasi-one-dimensional superconductor 



at zero magnetic field. Second, we calculate the quasi- 
particle density of states at zero temperature and the 
uniform spin susceptibility tensor at low temperatures 
for various candidate symmetries of the order parameter 
consistent with our group theoretical analysis. Third, 
we make connections to scanning tunneling microscopy 
(STM) of quasiparticle density of states and Knight shift 
measurements of the spin susceptibility tensor. 

We study quasi-one-dimensional systems with a single 
band, in an orthorhombic lattice, and allow for singlet or 
triplet pairing. We consider the following dispersion 



\t^\ cos(fc^a) - |ty| cos{kyb) - li^l cos(fc^c), (1) 



where \tx\ ^ \ty\ ^ \tz\. In the limit of weak in- 
teractions and low densities these quasi-one-dimensional 
systems exhibit a well defined Fermi surface which is 
open, being formed of two separate sheets which inter- 
sect the Brillouin zone boundaries. We work with the 
Hamiltonian H = H^in + Hint, where the kinetic en- 
ergy and the chemical potential contribution are H^in = 
I]k.a('^k - A')V'k.aV'k,a, and the interaction is 



Hint = I X! X! ^"/37<5(k,k')6i^(k,q)6^5(k',q) 



(2) 



kk'q a/37(5 



with &l^(k,q) = V'^k+q/2,aV^I+q/2,/3' "^^^^^ ^' ^ 

5 are spin indices and k, k' and q represent linear mo- 
menta. In the case of weak spin-orbit coupling and triplet 
pairing, the model interaction tensor can be chosen to be 

F„0^5(k,k') = r,^^5-l/r(k,k')<^r(k)0Mk'), 

where Vap-^s = v^p • with Vap = (icrcTy )q/3- In ad- 

dition, the interaction Vr corresponds to the irreducible 
representation F with basis function c/fr (k) representative 
of the orthorombic group. In the case of strong spin-orbit 
coupling the interaction 

Vc,p^s{k, k') = Vr{k, k') [$r(k) • v^^j] [*Mk') ■ v^5/2] 

where the interaction Vr corresponds to the irreducible 
representation F with basis function vector $r (k) repre- 
sentative of the orthorombic group. 
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In both weak and strong spin-orbit coupling, we can 
use either the equation of motion method or the 
functional integration method ||l^;|l3| in the zero center 
of mass momentum pairing approximation (which corre- 
sponds to the BCS limit in weak coupling) to obtain the 
anomalous Green's function 



Fa/3(k, iuJn 



Aa/3(k) 



and the single particle Green's function 



El 



(3) 



(4) 



where 



Tr 



ek — /i, /i is the chemical potential, E'k — 
is the quasiparticle excitation energy, and 

At(k)A(k) /2. The matrix A(k) has matrix 

elements AQ,^(k). The expressions for the single parti- 
cle (Eq. U) and for the anomalous (Eq. ||) Green's func- 
tions are valid only in the unitary case where A^(k)A(k) 
is diagonal. We will not discuss here the non-unitary 
case. Using the single particle and anomalous Green's 
functions defined above and standard many body meth- 
ods [^-|l3| we obtain the familiar forms 



A„,(k) Vf3.,5{K k')^^ tanh 



N 



(5) 



(6) 



where nqp(k) = (1 + ^kZ-^k) /(i?k) and n,/i(k) = 
(1 — ^k/-Ek) (1 — /(-Ek)) , for the order parameter and 
number equations, respectively. 

Next we consider the allowed symmetries of the order 
parameter Aq,^ for an orthorhombic crystal p^ , with a 
conventional symmetry normal state. Here, the relevant 
crystallographic point group is D2h, which has only one 
dimensional representations p5| ]. This implies that the 
order parameter matrix in the triplet channel. 



A(k) = i(At^(k)d(k) - a) a. 



J" 



(7) 



must transform according to the one dimensional repre- 
sentations of the orthorhombic point group D2h, under 
the assumption that the order parameter does not break 
the crystal translational symmetry, i.e., the order param- 
eter is invariant under all primitive lattice translations. 
Under the transformation k — k the three-dimensional 
vector d(k) is antisymmetric (odd), while the function 
Atr(k) is symmetric (even). 

Here, we will be interested in triplet states which do 
not break time reversal symmetry, and our analysis will 
be confined to zero magnetic field only. In Tables | and || 
we summarize the group theoretical analysis for A(k), 



in the weak and strong spin-orbit coupling cases respec- 
tively. The tables include the state nomenclature, the 
vector d(k), and the type of zeros of the quasiparticle ex- 
citation spectrum E^, when fl = fi — min[ek] is positive. 
In Table ^, the vector (0,0,1) is indicated up to an arbi- 
trary rotation in spin space. In Table ||, the numerical 
coefficients A, B and C are determined through Eqs. (|^) 
and (^. It is crucial to emphasize that the basis func- 
tions X(k),y(k) and Z(]i) transform like kx,ky and kz 
under the crystallographic point group operations. How- 
ever, these functions cannot be chosen to be equal to 
kx , ky and kz as done in the work by Lebed, Machida, 
and Ozaki ||l6| (LMO), since the Fermi surface e(k) = fi 
intersects the Brillouin zone along the y and z directions. 
Thus, it is necessary to take into account the periodic- 
ity of the order parameter matrix A(k) and of the order 
parameter vector d(k) in reciprocal (momentum) space. 
As a result, the minimal basis set must be periodic and 
may be chosen to be X(k) = sin(/cj,a), F(k) = sin{kyb) 
and Z{'k) = sin(fczc). 

For weak spin-orbit coupling (Table ||) the only candi- 
date for weak attractive interaction {p, > 0) is the state 
•^i?3„(a), where the quasiparticle excitation spectrum £'k 
has no zeros and is fully gapped. For strong spin-orbit 
coupling (Table |l|) there are three candidates for weaker 
attractive interaction (fl > 0), i.e., the states Aiu,Biu 
and B2u , where the quasiparticle excitation spectrum 
may have no zeros and may be fully gapped. When fl > 0, 
the state Aiu is fully gapped only for A and for any 
value of B and C; the state Biu is fully gapped only for 
B ^ and for any value of A and C; and the state B2u 
is fully gapped only for C and for any value of A 
and B. Note that in the case of strong attractive inter- 
actions where /i < 0, the excitation spectrum i?k is fully 
gapped jl^ for all states in both tables. 

TABLE I. Weak spin-orbit coupling. 



State 


d(k) 


£k = > 0) 




{0,0,1)XYZ 


lines 




(0,0, 1)Z 


lines 




(0,0,1)Y 


lines 




(0,0,1)^ 


none 


TABLE IL Strong spin-orbit coupling. 


State 


d(k) 


£k = (/i > 0) 


Am 


{AX, BY, CZ) 


none, points or lines 


Biu 


{AY, BX, CXYZ) 


none or lines 


B2u 


{AZ, BXYZ, CX) 


none or lines 


Bzu 


{AXYZ, BZ, CY) 


points or lines 
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Next we turn our attention to the calculation of 
the quasiparticle density of states (QDOS) and uniform 
spin susceptibility at low temperatures for two poten- 
tial triplet candidate states, which are fully gapped: a) 
the weak spin-orbit coupling state ■^i33„(a); b) the strong 
spin-orbit coupling states Aiu- jl^ 

The QDOS for these different symmetries can be ob- 
tained from the single particle Green's function as 



Af{u!) = Tr ImGap (k, iojn = uj + iS), 



(8) 



where GapiKii^n) is defined in Eq. (g). The QDOS, 
shown in Fig. 1 can be measured in STM experiments. 
Although these experiments have not yet been performed 
in (TMTSF)2C104 and (TMTSF)2PF6, our theoretical 
predictions can serve as qualitative guides for the extrac- 
tion of gaps and symmetry dependent features when ex- 
perimental results become available. In particular, STM 
measured gaps could be compared with gaps measured 
either thermodynamically (e.g., specific heat) or in trans- 
port experiments (e.g., thermal conductivity H). 






(d) 



12-2-1012 

co/Aj 

FIG. 1. Plot of QDOS versus frequency: the weak 
spin-orbit coupling state ^-B3„(a) is shown in (a); the strong 
spin-orbit coupling state Aiu is shown in (b) for A — 0.01, 
B = V2 A'2, C = 0; in (c) for A = B = 1, C = 0; in (d) 
for A = V2 - B"^, B = 0.01,C = 0. The parameters used are 



ioooa:, \u 



imK and 



5K, Atr = Ao = 3.0K 



and n = -250K. 



We compare in Fig. 1 the QDOS for the states 
^B^uic) (weak spin-orbit coupling) and Aiu (strong spin- 
orbit coupling) for various values of the constants A, 
B, and C. The symmetry dependent features of the 
QDOS are manifested through the magnitude of the 
order parameter vector d(k). For the ■^i33„(a) state 
|d(k)| cx |sin(fc2.a)|, while for the Ai„ state |d(k)| cx 
y/A^\sm{k^a)\^ + B^\ sin(fcy5)|2 + C^\ sin(fc^c)|2. In 
Fig. 1 (b), (c) and (d), we study only the case corre- 
sponding to C = 0, where A > _B, A = S, and A <^ B, 
respectively. Notice that Fig. 1(d) is nearly identical 
to Fig. 1(a) given that |d(k)| is essentially the same in 



this case. Furthermore, notice that while the position 
of the peaks in Figs. 1 (a), (b), (c) and (d) are essen- 
tially the same in scaled units (wp/Ao ~ ±1.40), the 
corresponding gap sizes in scaled units are respectively 
LUg/Ao — 1.32; 0.01; 0.94, 1.32. Gap sizes, peaks and the 
general shape of the QDOS should be in principle iden- 
tifiable in an STM experiment. However, such experi- 
ments alone cannot uniquely determine the symmetry of 
the order parameter in triplet quasi-one-dimensional su- 
perconductors, since the QDOS depends only on |d(k)|. 

Now, we turn our attention to the calculation of the 
spin susceptibility tensor, which explicitly depends on 
both the magnitude and direction of d(k), and, thus, 
may help ellucidate the symmetry of the order parame- 
ter in quasi-one-dimensional superconductors. The spin 
susceptibility tensor for a triplet superconductor is 

where we use the Einstein summation convention, and 
the four- vector = (q, w). The tensor {Pmn)ai3-yS = 
io'm)a0io'n)'yS, coutaius Pauli spin matrices, and the ten- 
sors 



Aafijsiqtj.) 



Fl^{-k + q, -it/ + iuj)Fi3s{K i^), 



contain the Green's functions described in Eqs. (|^) 
and d), and (3 = l/keT. For w ^ and q ^ 0, 



Xm"(0, 0) = ^ [Xmn,l(k) + (k)], 
k 

where the k-dependent tensors have the forms 
Xmn,i(k) = X||(k)i?e d,*„(k)d„(k), 

Xm«,2(k) = Xi(k) [Smn - Re d*j(k)d„(k) ) , 



(9) 



with (i„(k) — (i„(k)/|(i„(k)|. The parallel component is 



while the perpendicular component is 



XX (k) = 2^*1 



da 



2Ev 



(1 - 2/(i?k)) 



(10) 



(11) 



This result is more general than the expression quoted 
in LMO |l^ in a couple of ways. First, the expression 
derived in Eqs. (||), (|l^), and (jl^) includes particle-hole 
symmetry effects. Second they are valid at finite T. They 
do not include, however, Fermi liquid corrections. 
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In Fig. 2, we show the theoretical uniform Xmn only for 
the ^Bsu and Aiu states Q, where triangles correspond 
to Xii circles to X22, and squares to X33- It is known ex- 
perimentally (Knight shift) ||^ that the spin susceptibility 
of (TMTSF)2PF6 for H 1| b' is very close to xn- Experi- 
ments for magnetic field along other directions have not 
yet been performed. Thus, for definiteness, we choose 
the unit vectors 3 (m = 3), 2 (m = 2) and 1 (m = 1) to 
point along the b', a and c direction, respectively. 

For the orthorhombic symmetry Xmn is diagonal, and 
is calculated here under the assumption of constant d(k), 
i.e., the direction of d(k) is not changed upon applica- 
tion of a small magnetic field. In the case of state Aiu 
(strong spin-orbit coupling) a small magnetic field cannot 
rotate d(k) which is pinned to a particular lattice direc- 
tion. Here, Xmn is still diagonal, however the diagonal 
components are not equal in general. Thus, the experi- 
mentally measured Xmn and the theoretically calculated 
Xmn (^t constant d(k)) should agree for small enough 
magnetic fields. However, in the case of state ^B^uia) 
(weak spin-orbit coupling) a small magnetic field can al- 
ways rotate d(k) to be perpendicular to H, and thus min- 
imize the magnetic free-energy Fmag = -HmXmnHn/2. 
In this case, Xmn ~ XN^mn, where xn is the normal state 
value, for any direction H. Thus, Xmn ^nd Xmn (shown 
in Fig. 2) are different. 




0.10 



FIG. 2. Plot of the theoretical uniform spin suscepti- 
bility tensor components xn (triangles); X22 (circles); X33 
(squares) at low temperatures. The weak spin-orbit coupling 
state ^Bsu is shown in (a); strong spin-orbit coupling state 
Aiu is shown in (b) for A = 0.01, B = V 2 - A^, C = 0; in 
(c) for A = B = 1, C = 0; in (d) for A = V2 - B^, B = 0.01, 
C = 0. The parameters used are \t^\ = lOOOA", \ty\ = IQQK 
and 1 = 5A", At^ = Ao = S.OA' and ^ = -250if. 



In conclusion, we have studied order parameter sym- 
metry features in the quasiparticle density of states 
and spin susceptibility tensor of orthorhombic quasi-one- 
dimensional superconductors. We studied both the weak 
and strong spin-orbit coupling cases from a group the- 
oretical point of view at zero magnetic field. Based on 



experimental evidence, we would like to suggest that the 
weak spin-orbit coupling state ^i33„(a) is the best candi- 
date for the order parameter symmetry for these systems 
since this state is: (1) fully gapped and consistent with 
thermal conductivity measurements (2) characterized 
by weak spin-orbit coupling and consistent with weak 
spin-orbit coupling fits of Ta{H) for (TMTSF)2PF6 (l|] 
at low magnetic fields; (3) consistent with no observable 
Knight shift when H || 6', and predicted to have no ob- 
servable Knight shift for any direction of H. 

We would like to thank the Georgia Institute of Tech- 
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